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Abstract 

We investigate thermal entanglement between two non-nearest-neighbor sites in ferromagnetic 

' r-| ' 
Q I Heisenberg chain and on fractal lattices by means of the decimation renormalization-group (RG) 

method. It is found that the entanglement decreases with increasing temperature and it disappears 

beyond a critical value Tc- Thermal entanglement at a certain temperature first increases with the 

■^ I increase of the anisotropy parameter A and then decreases sharply to zero when A is close to the 

isotropic point. We also show how the entanglement evolves as the size of the system L becomes 
— L ■ large via the RG method. As L increases, for the spin chain and Koch curve the entanglement 

Q . between two terminal spins is fragile and vanishes when L > 17, but for two kinds of diamond-type 

o ■ 

hierarchical (DH) lattices the entanglement is rather robust and can exist even when L becomes very 
K^ large. Our result indicates that the special fractal structure can affect the change of entanglement 

1^^ . with system size 

rn 
O 
(N 
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I. INTRODUCTION 



An essential difference between quantum and classical physics is the possible existence 
of nonlocal correlation in quantum system which is called the entanglement [l|. Recently, 
the quantum entanglement has been recognized as an crucial resource in various fields of 
quantum information such as quantum communication and computation J2|-l7|. Since the 
entanglement is fragile and sensitive to many environment factors, many efforts are devoted 
to studying stable entanglement for realistic system in finite temperature. Thus, thermal 
entanglement has naturally received much attention by its advantage of stability and requir- 
ing neither measurement nor controlled switching of interactions in the preparing process. 
For spin systems can be used for gate operation in quantum computer thermal entanglement 



on solid spin systems have been wide 



in the absence 
clusters 



12 



y studied, for example, the Heisenberg spin chain both 



8| and presence (9|, [lO| of an external magnetic field , spin rings ll| and spin 



13| . However, most of these works only focused on thermal entanglement be 
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m- 



This motivate us 



tween nearest, next-nearest or next-to- next-nearest neighbor spins 
to propose two questions: (i) Can thermal entanglement exist between distant non-nearest- 
neighbor sites in spin system? (ii) How does thermal entanglement evolve as system size 
grows? But it is very difficult to obtain exact results on entanglement in spin systems on 
arbitrary lattices especially fractal lattices, since this usually requires the expression of the 
partition function which is too complicated to solve when the system size becomes very 
large. 

In recent years, the entanglement at zero temperature in the large size system has been 
studied by adopting the renormalizat ion-group (RG) method. In 2002, A. Osterloh et al first 
introduced the density matrix renormalizat ion-group (DMRG) approach to study the entan- 
glement close to the quantum phase transition (QPT) 20|] and reveal a profound difference 



between classical correlations and the non-local quantum correlation. Further,by applying 
the quantum renormalization-group (QRG) approach, M. Kargarian et al investigated the 



entanglement in the anisotropic Heisenberg model |2ll . |22| and discussed the nonanalytic 
behaviors and the scaling close to the quantum critical point of the system. Recently We 
have calculated the block-block entanglement in the XY model without and with staggered 



Dzyaloshinskii-Moriya (DM) interaction by using this QRG method and have found t 



he DM 



interaction can enhance the entanglement and influence the QPT of the system 



23 



24|. 



Inspired by above idea, we apply the real-space renormalization-group (RSRG) approach 



which is developed in the Refs. 25l-l35| to study the thermal entanglement between two end 
sites in the spin chain, Koch curve and on the diamond-type hierarchical (DH) lattices and 
analyze the influence of the temperature, the anisotropy parameter and the system size on 
the entanglement. 

II. MODEL AND METHOD 

The effective Hamiltonian of the spin-1/2 anisotropic ferromagnetic Heisenberg spin chain 
with L sites is 

L 

-PH = Y,K [(1 - A) {atat^, + 0^0^^^) + af af_,,] , (1) 

i 

where a'^ (a = x,y,z) denote the Pauh operators at site i. K = (3 J = J/ksT, J > is 
the exchange coupling parameter, and ks is the Boltzmann constant. For simplicity, we 
assume that fc^ = 1 and J = 1. The sum is over all the nearest-neighbor spin pairs and 
A G (— C)0, 1] is the anisotropy parameter. For A = and A = 1, the isotropic Heisenberg 
{XXX) and Ising model are obtained, respectively. The state of the above system in thermal 
equilibrium can be described by the density operator p = Z~^ (e~^^), where Z =Tr(e~^^) 
is the partition function. 

The entanglement of two-qubit system in the thermal state p^^ can be calculated by the 
negativity 



1 
36j which is based on the partial transpose method ^]. The negativity N is 



deflned as 

iV(p^J = 2^|/i,|, (2) 

i 

where /i, is the negative eigenvalue of (^ , and T\ denotes the partial transpose with respect 
to the flrst subsystem. According this deflnition one can easily obtain thermal entanglement 
A^ {K\ A') of the two-spin chain with the Hamiltonian iJ(2 (-^') A') (such as Fig 1. (a) n = 
shown). However, when the size of system becomes large, the density matrix is difficult or 
impossible to gain. The entanglement of two terminal spins on this system can not be 
directly worked out. 

We apply the decimation RSRG method to solve the above problem. This decimation 
RSRG method 27|, |29|, |38| has proved to be successful in spin chain and especially the fractal 



lattices. For the spin chain, this decimation procedure is illustrated in Fig. 1 (a). Simply, 
the generator is taken out in the infinite system. The generator with the Hamiltonian 
Hi32 {K, A) is renormalized into the new two-site chain with the Hamiltonian H[2 {K', A') 
by integrating the internal site 3 with the partition function being preserved. This transfor- 
mation can be described as 

exp {H[2) = Trs exp {H132) ■ (3) 

We can obtain the recurrence relation between the original parameters {K, A) and the new 



parameters {K',A') by solving the trace Tr3 with the method developed in Refs 29|, |31 |. 
Combining this relation K' = g (K, A) , A' = h {K, A) and the negativity A^ (K', A') of 
two-site system, the entanglement between two terminal spins in three-qubit chain can be 
obtained as follow 

N (K', A') = Nig (K, A) , h (K, A)) . (4) 

The entanglement between two distant terminal spins in the Heisenberg chain can be 
calculated after many iterations of the recurrence relation. We take use of the same method 
to study thermal entanglement between two terminal sites in Koch curve and on the DH 
lattices as shown in Fig. 1 (b) (c) and (d). The analytical expression about the negativity 
in Eq. (^ is difficult to obtain, we will show some numerical results. 

III. HEISENBERG SPIN CHAIN 

We first study how the entanglement between terminal sites in spin chain with differ- 
ent number of sites L varies with temperature T at A = —0.2 (shown in Fig. 2). For 
different cases of L, the results have the similar feature that the entanglement decreases 
monotonically with increasing temperature and it vanishes beyond the critical temperature 
Tc. At T = 0, the system is in the entangled ground state. As temperature increases, the 
entanglement decreases due to the mixture of the unentangled excited state with the ground 
state. At T = T^, the system is governed by the unentangled excited state completely, and 
therefore the entanglement vanishes. Comparing the entanglement of terminal sites in the 
different system, it is found that the entanglement of the ground state decreases sharply with 
increasing L. Different from the maximally entangled Bell ground state for L = 2 system, 
the ground state for L > 2 system becomes a degenerate and related to A state which cause 



the decrease of entanglement. The energy gap between the ground and the unentangled 
excited state increases with httle range but the thermal energy of system increases with 
large range when L increases. For L > 2, the system can easily overcome the gap and enter 
the unentangled state. That leads to the decrease of T^. This phenomenon reflects that 
thermal fluctuation of internal sites may suppress quantum effect. 

The influence of the anisotropic parameter A on the entanglement between two terminal 
sites at a fixed temperature T = 0.01 is plotted in Fig. 3. As can be seen, for L = 2 
case, the system firstly keeps in the ground state with the entanglement A^ = 1. Then, 
as A approaches the isotropic point A = 0, the energy gap between the ground and the 
unentangled excited states becomes so small that the system can jump to the unentangled 
excited state. Therefore there exists a sharp decease for entanglement and the entanglement 
vanishes when A is close to zero. This result also accords with that the entanglement can 



not exist in an isotropic Heisenberg ferromagnetic chain in Ref [39[. For L ^ 2 cases, it is 
found that the entanglement increases firstly with increasing A because the ground state is 
related to A and will change with A. All entanglement jump down to zero when A reaches 
to zero. From above results, we can see that the entanglement is fragile and when L > 17 
the entanglement does not exist whatever the temperature and the anisotropic parameter 
are. 

IV. FRACTAL LATTICES 

The properties of phase transition on different fractal lattices have been studied by the 
RSRG method, and the entanglement on these self-similar lattices remains to be explored. 
We turn to the study of the entanglement between end sites in Koch curves with non-integer 
fractal dimension df = In 4/ In 3 and plot the numerical results of negativity versus T and A 
for different L in Fig. 4. Compared with the entanglement in the spin chain (L = 5), it has 
similar properties that the entanglement decreases with T and the maximal of entanglement 
and Tc are approximately equal. But the entanglement variation versus A is very different 
from that in spin chain, i.e., the range that the entanglement can exist is smaller, the 
maximal of entanglement is lager. The entanglement decreases quickly as L becomes large 
and there exists no entanglement any longer when L > 17. We can deduce this result from 
the similar Hamiltonian and open boundary conditions of these two systems. 



Now we consider two kinds of DH lattices with fractal dimensions df = 2 (lattice A, 
for simplicity) and df = In 5/ In 2 (lattice B, for simplicity). The RG transformation on 
these DH lattices respectively have been shown in Fig. 1 (c) and Fig. 1 (d). We first 
discuss the dependence of the entanglement between terminal sites on T with A = —0.2. In 
Fig. 5, one can find some similar behaviors that the entanglement is the maximal value at 
T = 0, and the entanglement decreases with increasing temperature and vanishes beyond 
the critical temperature. However, some different phenomenons are also observed that the 
entanglement between end sites on these two lattices decrease more slowly with increasing 
L and it still exists even though the system size becomes very large {L = 1564). For the 
case of on lattice B, as L increases, the entanglement at zero temperature decreases but 
the corresponding Tc increases. It is obvious that the entanglement for different L crosses 
at T ?^ 0.49. This result indicates that these two fractal lattices have special energy level 
structure. The energy gap between the entangled ground state and the unentangled excited 
state is so large that the system can jump to the unentangled excited state only at higher 
T. 

The variation of the entanglement between end sites on DH lattices versus A at T = 0.01 
is also discussed. As can be seen in Fig. 6 (a) for the lattice A, the entanglement firstly 
increases as A increases, and then it quickly decays to zero when A reaches the isotropic 
point. The entanglement decrease very slowly when L becomes very large and there exist 
no cross point when A reaches zero. For the case of lattice B, Fig. 6 (b) shows that the 
entanglement also exhibits stable and it changes very little when A is not very close to zero. 
The entanglement mainly remains robust with the increase of L. In this graph, we also 
observe that a "entanglement crossing" occurs at A ?^ —0.045. At a fixed temperature, the 
thermal excited energy of the system is determined. Only when A is very close to zero, 
the energy gap between the entangled ground state and the unentangled excited state can 
become so small that the system can enters the unentangled state. It also indicates that 
the different fractal structure can influence the entanglement by changing the energy level 
structure of system. 



V. CONCLUSION 

We have investigated thermal entanglement between two end spins in Heisenberg chain, 
Koch curve and on two kinds of DH lattices with df = 2 and df = 2.32 by the decimation 
RG method. The effect of the temperature and the anisotropy parameter on thermal entan- 
glement is discussed. It is found that the symmetry of system and the thermal fluctuation 
can suppress or promote the quantum effect at different conditions. We also have noticed 
that the entanglement on some special lattices may exhibit different property when the sys- 
tem size L becomes large. The entanglement on two kinds of DH lattices is quite robust 
and it can survive even though L becomes very large in contrast to that in spin chain. The 
phenomenon of the "entanglement crossing" indicates that the special fractal structure does 
influence on the entanglement. 

Acknowledgments 

This work is supported by the National Natural Science foundation of China under Grant 
No. 10775088, the Shandong Natural Science foundation under Grant No. Y2006A05, and 
the Science foundation of Qufu Normal University. The authors would like to thank Hong- 
Bing Li, Yin-Fang Li, and Cong-Fei Du for fruitful discussions and useful comments. 



[1] 
[2] 
[3] 

[4] 
[5] 
[6] 

[7] 

[8] 

[9] 

[10] 



J. S. Bell, Physics (Lang Island city), 1, 195 (1964). 

Artur K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 

C.H. Bennett and S.J. Wiesner, Phys. Rev. Lett. 69, 2881 (1992). 

C.H. Bennett et al, Phys. Rev. Lett. 70, 1895 (1993). 

B. E. Kane, Nature (Londond) 393, 133 (1998). 

C. H. Bennett and D. P. Divincenzo, Nature (Londond) 404, 247 (2000d). 

M. A. Nielson and I. L. Chuang, Quantum Computation and Quantum, Inform,ation (Cam- 
bridge University Press, Cambridge, U.K., 2000). 
W. K. Wootters, Phys. Rev. A 63, 052302 (2001). 

M. C. Arnesen, S. Bose, V. Vedral, and S. Bose, Phys. Rev. Lett. 87, 017901 (2001). 
X. Wang, Phys. Rev. A 64, 012313 (2001). 



[11] X. Wang, Phys. Rev. A 66, 034302 (2002). 

[12] Indrani Bose and Amit Tribedi, Phys. Rev. A 72, 022314 (2005). 

[13] A. M. Souza, M. S. Reis, D. O. Soares-Pinto, I. S. Oliveira and R. S. Sarthour, Phys. Rev. B 

77, 104402 (2008). 
[14] K. Audenaert, J. Eisert, and M. B. Plenio R. F. Werner Phys. Rev. A 66, 042327 (2002). 
[15] L. Zhou, H. S. Song, Y. Q. Guo, and C. Li, Phys. Rev. A 68, 024301 (2003). 
[16] Guo-Feng Zhang, and Shu-Shen Li, Phys. Rev. A 72, 034302 (2005). 
[17] Zhao-Yu Sun, Kai-Lun Yao, Wei Yao, De-Hua Zhang, and Zu-Li Liu, Phys. Rev. B 77, 014416 

(2008). 
[18] Fardin Kheirandish, S. Javad Akhtarshenas, and Hamidreza Mohammadi Phys. Rev. A 77, 

042309 (2008). 
[19] A. M. Souza, M. S. Reis, D. O. Soares-Pinto, I. S. Ohveira, and R. S. Sarthour, Phys. Rev. B 

77, 104402 (2008). 
[20] A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature 416, 608 (2002). 
[21] M. Kargarian, R. Jafari, and A. Langari, Phys. Rev. A 76, 060304 (2007). 
[22] M. Kargarian, R. Jafari and A. Langari, Phys. Rev. A 77, 032346 (2008) 
[23] Fu-Wu Ma, Sheng-Xin Liu and Xiang-Mu Kong, Phys. Rev. A 83, 062309 (2011) 
[24] Fu-Wu Ma, Sheng-Xin Liu and Xiang-Mu Kong, Phys. Rev. A 84, 042302 (2011). 
[25] A. P. Young and R. B. Stinchcombe, J. Phys. C: SoUd State Phys 9, 4419 (1976). 
[26] R. Julhen, J. Fields and S. Doniach, Phys. Rev. Lett. 38, 1500 (1977). 
[27] M. Suzuki and H. Takano, Phys. Lett. A 69, 426 (1979). 
[28] H. Takano and M. Suzuki, J. Stat. Phys. 26, 635 (1981). 

[29] Anibal O. Caride, Constantino TsaUis, and SusanaL Zanette, Phys. Rev. Lett. 51, 145 (1983). 
[30] A. M. Mariz, R. M. Zorzenon Dos Santos, C. Tsalhs and R. R. Dos Santos, Phys. Lett. A 

108, 95 (1985). 
[31] A. M. Mariz, C. Tsalhs and A. O. Caride, J. Phys. C 18 4198 (1985). 
[32] S. R. White and R. M. Noack, Phys. Rev. Lett. 68, 3487 (1992). 
[33] J. Ricardo de Sousa, Phys. Lett. A 216, 321 (1996). 
[34] N. S. Branco, J. Ricardo de Sousa, Phys. Rev. B 62, 5742 (2000). 
[35] J. Ricardo de Sousa, N. S. Branco, B. Boechat and Claudette Cordeiro, Physica A 328, 167 

(2003). 



[36] G. Vidal, and R. F. Werner, Phys. Rev. A 65, 032314 (2002). 

[37] A. Peres, Phys. Rev. Lett. 77, 1413 (1996); M. Horodecki, P. Horodecki and R. Horodecki, 

Phys. Lett. A 223, 1 (1996). 
[38] Andre M. C. de Souza, Phys. Rev. B 48, 3744 (1993). 



[39] M. A. Nielsen, Ph. D. thesis. University of Mexico, 1198, e-print arXiv: quant-ph/0011036 



Figure captions: 

Fig. 1. The procedure of the RG transformation. From (a) to (d), it shows the trans- 
formation of one-dimensional spin chain, Koch curve, the diamond-type hierarchical lattice 
with fractal dimension df = 2 and df = In 5/ In 2. 

Fig. 2. The entanglement between two end sites in Heisenberg ferromagnetic chain versus 
temperature at A = —0.2 for different number of sites L (from top to bottom, L = 2, 3, 5 
and 9). 

Fig. 3. The entanglement between two end sites in Heisenberg ferromagnetic chain as 
anisotropy parameter A at T = 0.01 for different number of sites L (from top to bottom, 
L = 2, 3, 5 and 9). 

Fig. 4. The variation of the entanglement between end sites in Koch curve: (a) the 
entanglement versus temperature T at A = —0.2. (b) the entanglement versus A at T = 
0.01. 

Fig. 5. Upper panel: the negativity of two end sites on the DH lattice with df = 2 
versus T at A = —0.2. Lower panel: the negativity of two end sites on the DH lattice with 
df = In 5/ In 2 versus T at A = —0.2. The negativity for different value of L has a cross 
point at T ^ 0.49. 

Fig. 6. Upper panel: the evolution of the entanglement between terminal sites on DH 
lattice with df = 2 as A increases at T = 0.01 for different value of L. Lower panel: the 
entanglement between terminal sites on DH lattice with df = In 5/ In 2 versus A at T = 0.01 
for different value of L. The entanglement for different L has a cross point at A ?^ —0.045. 
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